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Abstract. We study the Wigner functions of the nucleon which provide muhidimensional images of 
the quark distributions in phase space and combine in a single picture all the information contained 
in the generalized parton distributions (GPDs) and the transverse-momentum dependent parton 
("^ ' distributions (TMDs). In particular, we present results for the distribution of unpolarized quark 

^\i I in a longitudinally polarized nucleon. Treating this distribution as it was a classical distribution, 

we also obtain the expectation value of the orbital angular momentum operator and compare 
the corresponding results from different quark models with the values obtained using alternative 
definitions of the orbital angular momentum, as given in terms of the GPDs and the TMDs. 
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INTRODUCTION 



■^ . Quark Wigner distributions open a new way to map the distributions of momentum and 

Q^ I spin of the proton onto its constituents. They provide joint position-and-momentum (or 

cn ' phase-space) distributions, encoding in a unified picture the information obtained from 

(^ . transverse-momentum dependent parton distributions (TMDs) and generalized parton 

O I distributions (GPDs) in impact-parameter space. The concept of Wigner distributions 

in QCD for quarks and gluons was first explored in Refs. [1,2]. Neglecting relativistic 
effects, the authors introduced six-dimensional (three position and three momentum co- 
ordinates) Wigner distributions. In a recent work [3], we used the connection between 
rS I Wigner distributions and generalized transverse-momentum dependent parton distribu- 

te ■ tions (GTMDs) [4] to study five-dimensional distributions (two position and three mo- 

mentum coordinates) as seen from the infinite-momentum frame (IMF). The main ad- 
vantage of working in the IMF is that one obtains Wigner distributions which are com- 
pletely consistent with special relativity. However, it is well known that the phase- space 
distributions do not have a density interpretation because the uncertainty principle pre- 
vents to determine simultaneously position and momentum of a quantum-mechanical 
system. Accordingly, Wigner distributions are not positively defined. Nevertheless, the 
physics of the Wigner distributions is very rich and one can select certain situations 
where a semiclassical interpretation is still possible. 

The purpose of this contribution is to investigate the phenomenology of the quark 
Wigner distributions. As a matter of fact, since it is not known how to access these 
distributions directly from experiments, phenomenological models are very powerful in 
this context. Collecting the information that one can learn from quark models which 



were built up on the basis of available experimental information on GPDs and TMDs, 
one can hope to reconstruct a faithful description of the physics of the Wigner distribu- 
tions. To this aim we will rely on models for the light-cone wave functions which have 
already been used for the description of the GPDs [5, 6], the TMDs [6, 7, 8, 9, 10] and 
electroweak properties of the nucleon [6, 11, 12, 13, 14]. 

WIGNER OPERATORS AND WIGNER DISTRIBUTIONS 

Similarly to Refs. [1, 2], we define the Hermitian Wigner operators for quarks at a fixed 
light-cone time );+ = as follows 

with y^ = [0,0,b±], p+ the average nucleon longitudinal momentum and x = k'^ /p^ 
the average fraction of nucleon longitudinal momentum carried by the active quark. The 
superscript F stands for any twist-two Dirac operator F = 7+, 7^75, /cr-'+ys with j = 1,2. 
A Wilson line W = W {y — j,y + j\n) ensures the color gauge invariance of the Wigner 
operator. In the following we will focus on the quark contribution, ignoring the gauge- 
field degrees of freedom and therefore reducing the gauge link W to the identity. 

We define the Wigner distributions p^^\b^^,k^^,x,S) in terms of the matrix elements 
of the Wigner operators (1) sandwiched between nucleon states with polarization S. As 
outlined in Ref. [3], such matrix elements can easily be interpreted as two-dimensional 
Fourier transforms of the GTMDs in the impact parameter space. Although the GT- 
MDs are in general complex- valued functions, their two-dimensional Fourier transforms 
are always real-valued functions, in accordance with their interpretation as phase-space 
distributions. We note that bj^ and k± are not Fourier conjugate variables, like in the 
usual quantum-mechanical Wigner distributions. However, they are subjected to Heisen- 
berg's uncertainty principle because the corresponding quantum-mechanical operators 

do not commute [Sj_,^i_] 7^ 0. As a consequence, the Wigner functions can not have 
a strict probabilistic interpretation. There are in total 16 Wigner functions at twist-two 
level, corresponding to all the 16 possible configurations of nucleon and quark polar- 
izations. Here we will discuss only one particular case, namely the distortion of the 
distribution of unpolarized quarks due to the longitudinal polarization of the nucleon 

PLU = P^^^\i>_i,kj^,x,+ez) - p^^'^\b^^,kj_^,x,-ez). 

In Fig. 1, the upper panels show the distortions in impact-parameter space for u 

(left panels) and d (right panels) quarks with fixed transverse momentum k± = k± Cy 
and k^ = 0.3 GeV, while the lower panels give the corresponding distortions in the 
transverse-momentum space with fixed impact parameter S^ = b^Cy and b^ = 0.4 fm. 

We observe a clear dipole structure in both these distributions, with opposite sign for 
u and d quarks. We learn from these figures that the orbital angular momentum (0AM) 
of u quarks tends be aligned with the nucleon spin, while the 0AM of d quarks tends be 
antialigned with the nucleon spin. In particular, we notice that the distortion induced by 
the quark 0AM is stronger in the central region of the phase space {k± <^ and b± <^), 
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FIGURE 1. The distortions of the transverse Wigner distributions of unpolarized quarks due to the spin 
of the proton (pointing out of the plane). Upper panels: distortions in impact-parameter space with fixed 
transverse momentum kx_ = k±ey and k± = 0.3 GeV. Lower panels: distortions in transverse-momentum 

space with fixed impact parameter b± = b±ey and b± — 0.4 fm. The left (right) panels show the results 
for u {d) quarks. 



for both u and d quarks. The distortion in the b± space (see upper panels of Fig. 1) is 
more extended for d quarks than for u quarks, whereas the opposite behavior is found 
for the distortion in the Ix space (see lower panels of Fig. 1). In the case of 6? quarks, we 
also observe a sign change of the distributions in the outer regions of phase space {k± ^ 
and b± ^) which corresponds to a flip of the local net quark 0AM. By treating the 
Wigner functions as if they were classical distributions, we can obtain the expectation 
value of the quark 0AM operator by calculating the integral over the phase space of the 
distribution in Fig. 1 multiplied by (b± x k_i_)z- This definition differs from the quark 
0AM calculated from GPDs, using the Ji's sum rule [15], as well as from the quark 
0AM defined in terms of the hjj TMD [16, 17]. However, it turns out that in models 
without gauge-field degrees of freedom all the three definitions give the same values for 
the total quark contributions, while they differ for the individual u and d contributions. 
In particular, we present in Table 1 the results from the LCCQM and the light-cone 
version of the chiral quark-soliton model (;^QSM) restricted to the three-quark sector 



[11, 12, 13]. In the LCCQM, there is more net quark 0AM (£^Lf = 0.126) than in the 

TABLE 1. The results for quark OAM obtained from the Wigner 
functions (£^), from the Ji's sum rule (Lj.), and from the hfj TMD 
(^/) within the LCCQM and the jQSM for u-, d- and total {u + d) 
quark contributions. 



Model 

q 



LCCQM 

d TOT 



ZQSM 

d TOT 



tl 


0.131 


-0.005 


0.126 


0.073 


-0.004 


0.069 


Li 


0.071 


0.055 


0.126 


-0.008 


0.077 


0.069 


^; 


0.169 


-0.042 


0.126 


0.093 


-0.023 


0.069 



;t:QSM (Y^qLl = 0.069). For the individual quark contributions, both the LCCQM and 
the ;t:QSM predict that £f (from the Wigner functions) and ^/ (from the hjj TMD) 
are positive for u quarks and negative for d quarks, with the w-quark contribution larger 
than the J-quark contribution in absolute value. For Lf (from the Ji's sum rule) the 
LCCQM predicts the same positive sign for the u and d contributions, with the isovector 
combination L" —Lf > 0, similarly to a variety of relativistic quark model calculations. 
Instead, the ;^QSM gives L" < and Lf > 0, and therefore L^ — Lf< 0, in agreement 
with lattice calculations [18]. 

In summary, we have presented the first model calculation of the Wigner distributions 
within the light-cone formalism, using a derivation which is not spoiled by relativistic 
corrections. The results within a light-cone constituent quark model and the light-cone 
version of the chiral quark-soliton model are very similar, and allowed us to sketch some 
general features about the behavior of the quarks in the nucleon when observed in the 
bx_ plane at fixed k^, or in the ^^ plane at fixed b^. 
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